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Abstract 

We investigate the effect of enhancing gravity on saturated nucleate pool boiling 
of oxygen for effective gravities of lg, 6-0g, and 16g (g = 9.8 ms -2 ) at a satura- 
tion pressure of 760 torr and for heat fluxes of 10 ~ 3000 Wm~ 2 . The effective 
gravity on the oxygen is increased by applying a magnetic body force generated 
by a superconducting solenoid. We measure the heater temperature (expressed as 
a reduced superheat) as a function of heat flux and fit this data to a piecewise 
power-law/linear boiling curve. At low heat flux (< 400 Wm~ 2 ) the superheat is 
proportional to the cube root of the heat flux. At higher heat fluxes, the superheat 
is a linear function of the heat flux. To within statistical uncertainties, which are 
limited by variations among experimental runs, we find no variation of the boiling 
curve over our applied gravity range. 
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1 Introduction 



Nucleate pool boiling is one of the enduring problems in fluid mechanics and 
thermal physics. The large number of physical variables, many of them not 
directly measurable or controllable, and the difficulty in reproducing exact 
experimental conditions have thwarted attempts at constructing an accurate 
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theoretical description, or even phenomenological description, of this everyday 
process. In an attempt to shed more light on this difficult topic, we have 
chosen to study the effect of varying one particular parameter which has been 
inaccessible to most investigators: gravity. 



Many widely-used empirical correlations of boiling heat transfer neglect grav- 
ity. Several make no mention of gravity at all (Cooper [lj], Stephan & Abdelsal- 
am[2|). Those that do include the gravitational acceleration g (Rohsenow[3t], 
Forster & Greiffl), often do so only as a dimensional constant, not as a phys- 
ical variable. In these models varying the value of g yields unphysical predic- 
tions, as was recognized by DhirpJ, who suggests that these correlations are 
only valid when g is held constant at its Earth value of 9.8 ms~ 2 , regardless 
of the actual gravity in the experiment. To develop adequate models which 
account for gravity variations, more experimental data are needed. 



Experiments at high effective gravities using centrifuges [a, |7| and parabolic- 
trajectory aircraft [8j show results which vary depending on the applied heat 
flux. At low heat fluxes, increasing gravity increases heat transfer. At high 
heat fluxes, however, increasing gravity decreases heat transfer 0. 



Low gravity studies show an even wider range of results. The general con- 
sensus among investigators is that the factors influencing nucleate boiling can 
be qualitatively divided into macroscopic and microscopic phenomena. Macro- 
scopic effects, such as buoyancy, bubble dynamics, and thermocapillary effects, 
tend to depend on gravity. Microscopic effects at the heater surface, such as 
intermolecular forces, microlayer evaporation and microlayer conduction, are 
independent of gravity. As gravity is reduced, the microscopic phenomena 
begin to dominate over the macroscopic effects and compete with one an- 
other, resulting in varying experimental results depending on the details of 
the experiments 0, HQ. 



In this article, we investigate the effect of enhanced gravity in our Earth- 
based laboratory; to do this we must apply some additional body force to 
the fluid. Previous studies have used mechanical forces to create this body 
forcejg 1 , Q, However, we chose to use the large magnetic susceptibility of 
oxygen to produce the additional body force by applying a relatively modest 
magnetic field. By utilizing existing equipment in our laboratory, in particular 
a 10-cm bore, 9-T solenoidal superconducting magnet, we have constructed an 
inexpensive test apparatus for investigating a range of net effective gravities 
g e s between microgravity and 90g, with g = 9.8 ms~ 2 representing the Earth's 
gravitational acceleration. Here we present data for saturated nucleate pool 
boiling under modestly increased gravities, between lg and 16 g of net effective 
gravity. 
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2 Basic principle 



We simulate enhanced gravity by applying a body force to the test fluid that 
is the result of the interaction of the fluid's magnetic susceptibility with an 
applied magnetic field. The magnetic acceleration a on a fluid with magnetic 
susceptibility \ in a magnetic field B can be expressed as 

a=--(B-V)B, (1) 

Po P 



where p is the mass density of the fluid, /io is the vacuum permeability, and all 
quantities are in SI units. For oxygen near the normal boiling point (90.2 K) 
the ratio x/p equals +3.02 x 10~ 6 m 3 kg _1 for both liquid and gas phases fl2|] . 
We assume any temperature dependence in \l P is negligible over our tem- 
perature range of 90 — 95 K. To achieve lg of acceleration on oxygen at the 
normal boiling point requires |(B • V)B| = 4.5 T 2 m _1 . The magnetic force 
causes an acceleration that is independent of the fluid density and therefore 
independent of the phase of the fluid. 

To produce our magnetic field we use a superconducting solenoid manufac- 
tured by American Magnetics, Inc. (Oak Ridge, TN, USA). The field of the 
magnet is approximately that of a thin uniformly-wound finite solenoid. In 
particular, the on-axis magnetic field (shown in Figure la) is given by 

^ + «/2 z-a/2 \ 

^{z + a/2Y + V yj(z - a/2) 2 + 6 2 ) ' 



with parameters a = 25.5 cm, b = 7.78 cm and C = 6.28 x 10~ 2 T A -1 (deter- 
mined by fitting to data provided by the manufacturer), / being the current 
in the magnet measured in amperes (maximum of 84 A for this magnet), and 
the axial position z measured relative to the center of the magnet. The corre- 
sponding on-axis applied acceleration on the oxygen is plotted in Figure lb. 
The off-axis fields at the heater surface [z = 8.8 cm), calculated using formu- 



las in Ref. [13j], and the resulting accelerations are plotted in Fig. lc and Id, 



respectively. We ignore the radial force in our analysis. 



To obtain an approximately flat force profile over the volume of the liquid, 
we place the sample cell near the end of the magnet where the force obtains 
its maximum magnitude, with the center of the sample cell 11.0 cm above 
the magnet center, as illustrated in Fig. 1. Future studies at reduced effective 
gravity may be performed by placing the sample at the bottom end of the 
magnet, where the magnetic force counteracts the Earth's gravity. 
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3 Apparatus 



A schematic of our experimental setup is shown in Figure 2a. The sample cell 
is a quartz cylinder (50 mm O.D., 2 mm thickness, and 43 mm height, inner 
volume of 70 cm 3 ) with a stainless steel top plate and an oxygen- free high 
conductivity (OFHC) copper bottom plate (area of 16.6 cm 2 ), joined together 
by indium gaskets held in compression with 6-32 bolts tightened to 5 N • m of 
torque. Lakeshore model CX-1050-SD Cernox temperature sensors are placed 
on the interior surface of the top plate, in the center of the cell, and embedded 
in the bottom plate. The bottom sensor measures the heater plate temperature 
and is the primary temperature sensor used for analysis of the experiment. 

The top and center temperature sensors are used only during sample filling. 
During the boiling phase of the experiment the temperatures measured by 
these sensors drift depending on the liquid height and therefore are not useful 
for quantitative analysis. 

For the heat source, a 75-fi Minco resistive foil heater is attached to the 
exterior of the copper plate using Apiezon N grease and a beryllium copper 
clamp. The interior surface of the copper plate is roughened using 200-grit 
sandpaper (to provide nucleation sites) and then cleaned in consecutive washes 
of detergent, acetone, and isopropanol. Between experiments the sample cell is 
kept under vacuum to prevent adsorption of gas onto the heater surface which 
could change the boiling dynamics. To measure the magnetic field, Toshiba 
THS-122 Hall-effect sensors are located on each end of the sample cell, on-axis. 
A fiberoptic borescope built by Myriad Fiber (Dudley, MA, USA) is positioned 
adjacent to the quartz sample cell wall to provide imaging of the experiment. 
The borescope image is output via a Panasonic GP-KS162HD CCD camera 
to a television monitor. The temporal and spatial resolution of the camera 
are insufficient for quantitative measurements of bubble dynamics but serve 
nicely for diagnostic purposes and qualitative observations such as the level of 
the liquid during filling and identification of bubbles to confirm boiling. 

The heater and the bottom temperature sensor are controlled by a Lakeshore 
LS-340 temperature controller. The other temperature sensors and the mag- 
netic field sensors are measured in a 4- lead configuration using Keithley 182 
voltmeters and Keithley 220 current sources. 

A 6-mm O.D. stainless steel pumpline extends from the top of the sample cell 
to room-temperature gas handling equipment. Pumping of the oxygen sample 
gas is performed by Varian Vacsorb sorption pumps, filled with zeolite get- 
ter and cooled by liquid nitrogen. In consideration of the safety challenges of 
working with oxygen gas at potentially high pressures, we chose these pumps 
for their low cost and ease of use relative to fluorinated lubricant filled me- 
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chanical pumps. An MKS 65 1C pressure regulator, actuating an MKS 253 A 
butterfly value, controls the system pressure, which is measured by an MKS 
626A Barotron sensor placed in front of the butterfly valve. 

All of the above electronic sensing and control equipment are connected via 
GPIB to a PC running LabVIEW 7.0 software. The sensors are polled se- 
quentially with a cycle of approximately 10 seconds. The control relationships 
between the various components are summarized in Figure 2b. 



4 Procedure 

We prepare an experimental run by first precooling the magnet dewar with 
liquid helium and precooling the sample probe in a separate liquid nitrogen 
dewar. Approximately 0.1 torr of helium gas is added to the vacuum space of 
the probe to facilitate heat exchange between the sample cell and the nitrogen 
bath during filling. We fill the sample cell by introducing room-temperature 
oxygen gas into the sample cell via the pumpline. The gas condenses until the 
cell is filled with liquid oxygen, a process which takes about 90 minutes. The 
liquid level can be monitored visually using the borescope and can also be 
inferred by monitoring the various temperature sensors inside the cell. After 
the oxygen fills the sample cell, we close the fill valve and remove the helium 
exchange gas from the vacuum space using a turbomechanical pump, typically 
pumping overnight to a pressure of 10~ 6 torr. The probe is then moved to the 
liquid-helium-filled magnet dewar. As the probe exterior cools to liquid helium 
temperature, a charcoal pot attached to the interior surface of the vacuum can 
absorbs any residual exchange gas in the vacuum space. 

We begin monitoring the temperature and magnetic field sensors at this time. 
Figure 3 contains heater power, heater temperature, and system pressure data 
taken during a typical experimental run. Details of the various stages of the 
experiment are described below. 

To prevent the sample from freezing during the initial setup of the experiment 
(region "a" of Fig. 3), we maintain the heater temperature (measured by the 
sensor embedded in the copper plate, see Fig. 2a) at 90 K (just below the 
1 atm boiling point) using a Lakeshore LS-340 temperature controller. We 
record the heater power required to maintain this constant temperature with 
all valves closed, so that the system volume is fixed, and interpret this value 
as heat lost to the liquid helium bath by conduction through the support 
structures of the probe and the residual helium gas in the vacuum space; we 
subtract this quiescent heater power (typically, ~ 50 mW) from later heat 
measurements. While the temperature of the sample probe is stabilizing, we 
increase the magnet current to the desired value. 
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To begin the experiment proper, we open the valves to the sample and set 
the pressure controller to the desired value (760 torr for the data presented 
here). We disable the temperature control loop and set the heater power to 
the desired maximum value for the run (5.0 W for the data presented here, 
corresponding to a maximum heat flux of 3.0 kW m~ 2 ). As the temperature of 
the sample increases, the transition from natural convection to boiling appears 
as a sudden drop in the heater temperature because of the increased heat 
transfer efficiency in the boiling state (see the dashed line at the beginning 
of region "b" in Fig. 3, although see the discussion below). The temperature 
then continues to increase until the system reaches some initial steady-state 
with approximately constant temperature and pressure. 

After the system stabilizes at this initial heater power, we decrement the 
power by steps, typically in geometric series (region "c" in Fig. 3). At each 
new value of the power, we wait for the system to stabilize (about 5 minutes) 
and record the new steady-state heater temperature. The experiment ends 
when the heater power decreases below the threshold to maintain a constant 
temperature. This threshold heater power typically agrees within 10% of the 
quiescent heater power measured during the setup phase of the experiment. 

The time evolution of temperature and pressure during the initial warm-up 
(region "b" in Fig. 3) varies between runs because of transient behavior in the 
pressure regulation control loop. A common occurrence is that as the system 
transitions from natural convection to boiling (characterized by a drop in 
heater temperature), the sudden increase in gas flow upsets the feedback loop 
of the pressure controller, causing large pressure and temperature oscillations 
for several minutes while the controller re-establishes steady-state behavior. 
Some consequences of this transient behavior are discussed in the sections 
below. 



5 Analysis and discussion 



To characterize our data, we first rescale the heater temperature T to a reduced 
superheat t relative to the saturation temperature of the oxygen Tj,: 

t = (3) 



The saturation temperature T b is calculated by measuring the system pressure 
on the external pum pline (see Fig. 2b) and then applying the equilibrium 
equation of state [141] . To directly measure the saturation temperature would 
require placing temperature sensors at the liquid-gas interface. Because the 
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liquid level changes throughout the experiment our fixed-position sensors are 
inadequate for this measurement. 



The purpose of changing variables to the reduced superheat t is to provide 
a dimensionless parameter for analysis and facilitate comparison with future 
measurements at varying saturation pressures. We anticipate, based on some 
preliminary measurements not presented here, that the boiling curves at dif- 
ferent saturation pressures will collapse onto a single curve when expressed in 
these units. 

The reduced superheat we measure is well-fit by a piecewise function of the 
applied heat flux q: 



where q t is some transition heat flux (found by fitting), and the dimensionless 
quantities a, b, m, and c are fitting parameters subject to smoothness con- 
straints that t(q) and dt(q)/dq. The constraints are equivalent to 



leaving three independent free parameters per fit. Experimental data sets and 
the resulting fits are shown in Figure 4 with the fitting parameters listed in 
Table 1. The power law form was chosen for comparison with many empirical 
and theoretical correlations in the literature [H, HI, 0, 0]. However, this form 
could not fit the full range of our data. We have appended the linear portion of 
the fitting function and smoothness constraints to remedy the poor fits at high 
heat flux without adding an excessive number of free parameters. Averages 
for each value of g e s are also shown in Table 1. Because of irreproducibility 
between runs, the uncertainties in the average values are much larger than the 
uncertainties of individual runs. The suspect this irreproducibility is caused 
by uncontrolled transient effects, as discussed below. 

The piecewise boiling curve represented by Eq. 4 agrees qualitatively with 
the results of the centrifuge experiments reported in Refs. ja, 0], showing dif- 
ferent behavior at low and high heat fluxes. In particular, our low flux data 
approximately obey the power law t oc q 1 ^ 3 to within statistical uncertainties; 
this form has precedent in some of the more well-known correlations, such as 



Cooper Stephan & Abdelsalam[2], and Rohsenow[3j. This sublinear behav- 
ior may indicate that the number of nucleation sites changes over this range of 
heat flux. On possibility is that q t represents the minimum heat flux required 




(4) 



m = ab; 
c— a(l — b), 



(5) 
(6) 
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Fig. 1. Plots of the magnetic field and resulting acceleration for a magnet current 
of 20 A. Panel a: On-axis magnetic field B z as a function of axial position z. Panel 
b: On-axis magnetic acceleration a z on oxygen as a function of axial position z. 
Panel c: axial component of magnetic field B z (dashed line) and radial component 
B r (solid line) evaluated at the heater surface (z = 8.8 cm) as a function of distance 
from the axis r. Panel d: axial component of magnetic acceleration a z (dashed line) 
and radial component a r (solid line) evaluated at the heater surface (z = 8.8 cm) 
as a function of distance from the axis r. The additional —lg of acceleration in the 
— z direction from the Earth's gravity is not shown in these plots. The approximate 
location of the sample cell (4.4 cm in height, centered at z = 11.0 cm, and 2.3 cm 
in radius) is indicated by the shaded areas. Not shown is the size of the magnet, 
which extends from —15 cm < z < +15 cm. 



Run # 


9es/g 


?(/(Wm- 2 ) 


a/10- 3 


b 


m/10~ 3 


c/10~ 3 


1 


1 


1135 ± 12 


6.58 ± 0.05 


0.572 ± 0.009 


3.77 ±0.03 


2.81 ±0.06 


2 


1 


729 ±6 


7.96 ± 0.03 


0.363 ± 0.005 


2.89 ±0.02 


5.07 ±0.03 


Average 


1 


810 ± 160 


7.59 ± 0.66 


0.412 ± 0.090 


3.16 ±0.41 


4.62 ± 0.89 


3 


6.0 


272.7 ±0.7 


6.27 ±0.01 


0.2586 ± 0.0008 


1.621 ±0.003 


4.65 ±0.01 


4 


6.0 


497.4 ± 1.1 


6.55 ±0.01 


0.3983 ±0.0012 


2.609 ± 0.004 


3.94 ±0.01 


Average 


6.0 


340 ± 100 


6.41 ±0.14 


0.302 ±0.064 


1.98 ±0.47 


4.30 ±0.29 


5 


16 


586 ±2 


9.28 ±0.02 


0.386 ±0.001 


3.585 ± 0.006 


5.69 ±0.01 


6 


16 


940 ±4 


11.27 ±0.02 


0.317 ±0.001 


3.574 ± 0.009 


7.69 ±0.01 


Average 


16 


660 ± 140 


10.28 ±0.94 


0.352 ±0.029 


3.582 ±0.005 


6.7 ±1.0 


Overall 


average 


370 ± 150 


7.2 ±1.6 


0.33 ±0.06 


2.3 ±0.8 


5.5 ± 1.2 



Table 1 

Fitting parameters for the data in Fig. 4, showing effective gravity <7 e ff , and param- 
eters from Eq. 4. Average values for each effective gravity and for the entire data 
set are also shown. 
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to maintain activity each individual nucleation site. As the heat flux decreases 
below this threshold, these sites begin to close off, resulting in decreased heat 
transfer efficiency. This exponent does not appear to vary with g e g, at least 
to within the uncertainties of our data. 

We are unclear about what causes the t oc q behavior in our high flux data, 
although the shape of this part of the curve is unmistakable in the linear-scale 
data plots (see the left panel of Fig. 4). Within the simple picture of nucleation 
sites discussed the previous paragraph, this is consistent with the conjecture 
that the number of nucleation sites is independent of heat flux above the 
threshold flux q t . In the energy added by increasing heat flux goes only into 
increasing the superheat, rather than into creating new nucleation sites. 

The functional form we have chosen clearly appears in all of our data sets, 
and the fits for each individual data set are good with the uncertainties in 
the fitting parameters being a few parts per thousand. However, because of 
irreproducibility the averages of runs have large uncertainty. We have taken 
some steps to minimize this irreproducibility. 

One technique which improves reproducibility is to avoid starting the sys- 
tem in the metastable natural convective state by starting at a maximum 
heater power and then decreasing the power. This sequence minimizes any ir- 
reproducibility caused by the hysteretic nature of the transition from natural 
convection to boiling. The system seems to retain some memory of the maxi- 
mum heater power applied over any given experimental run. To demonstrate 
this memory effect we produced a pair of runs (Fig. 5) in which the heater 
power was alternated. In the first data set, shown in Fig. 5a, the heater power 
starts at a nominal maximum value, steps to a lower value, returns to the 
maximum value, steps down to a still lower value, and the pattern repeats. 
The key effect in this data is that the heater temperature returns to the same 
value with each return of the power to the maximum power, indicated by the 
dashed line in the lower panel of Fig. 5a. Also, the sequence of lower power 
steps result in a monotonically decreasing heater temperature sequence. This 
pattern is in contrast to what happens if the heater power starts at some low 
value, increases, returns to the low value, and then increases to a yet higher 
value. Such a data set is shown in Fig. 5b. The important effect in this data is 
that the temperature at the 2 W steps (indicated by the arrows in the figure) 
seems to depend on the previous heat value, unlike Fig. 5a. 

We hypothesize that the number of nucleation sites on the heater surface is de- 
pendent on the maximum heater power applied since the initiation of boiling. 
This could be verified by imaging the heater surface and counting nucleation 
sites, but our current apparatus lacks this capability. To ensure the minimum 
variation between experimental runs, we've chosen a maximum nominal heater 
power of 5 W, corresponding to a maximum heat flux of 3 kW m -2 . At heater 
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powers higher than 5 W the liquid boils off too quickly for us to acquire suf- 
ficient data, exhausting our small sample volume before the system can reach 
a steady state. 

Unfortunately, some irreproducibility between runs persists, most likely be- 
cause of the transient fluctuations during the preparation phase of the indi- 
vidual data runs (for example, region "b" in Fig. 3). We believe these oscilla- 
tions (or occasional lack thereof) cause the remaining variability of the data. 
In particular, we observe that the steady-state temperature of the heater at 
the maximum heater power varies between runs (the rightmost points of the 
boiling curves in Fig. 4). Unfortunately, in some runs the transient pressure 
drops below the minimum range of our pressure gauge (0.5 torr), so we are un- 
able to quantify any possible relation between the minimum transient pressure 
and the heat transport data. We suspect the variation in maximum temper- 
ature between runs is evidence of a variation in the number of nucleation 
sites. When the pressure drops during these transient oscillations, the satu- 
ration temperature decreases causing a temporary superheated condition in 
the liquid; the sudden excess of energy in the system may go into activating 
nucleation sites. These new nucleation sites may remain active after the sys- 
tem reaches a steady state, resulting in more efficient heat transfer (therefore, 
lower heater temperature) than would result if the pressure behaved smoothly 
during the initial stages of boiling. 

The irreproducibility caused by the above transient effects may be also obscur- 
ing any gravity effect on heat transport. The values of the fitting parameters 
given in Table 1 are independent of the applied force to within statistical un- 
certainties. This is consistent with previous experiments with gravity values 
in our range 1 < g e s/g ^ 16 which also find boiling heat transfer to be in- 
dependent of gravity jg, Q, @|. The large uncertainties in the average fitting 
parameters are the results of variations among data sets, corresponding to 
the difficulty in controlling the transient behavior during the setup of the ex- 
periment. A better understanding of the transient pressure effects and/or a 
larger number of data sets are probably required to resolve any gravity effect 
in the data. Comparison of parameter values among the various runs in Table 
1 also suggests that some parameters may be correlated, for example, qt and 
a, although more investigation is needed to confirm any such relation. 



6 Conclusions 

For the saturated nucleate boiling of oxygen, we observe two heat transport be- 
haviors. At high heat flux values q > 400 Win" 2 , the superheat varies linearly 
with the applied heat flux. At lower heat fluxes, the superheat is proportional 
to g 1 / 3 . The latter result is consistent with well established correlations from 
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the literature [H, 0, 0]. We are unable to see any change in heat transport 
behavior caused by varying effective gravity over the range lg < g c & < 16g, 
although our analysis is limited by some residual irreproducibility between ex- 
perimental runs related to transient pressure behavior during the setup phase 
of the experiment. 

Our apparatus has the capability to extend the range of gravity values studied 
t° < 9es/g ^ 90. The larger values of effective gravity can be obtained by 
increasing the current in our magnet up to its maximum value of 84 A. To 
study values of g e ff less than Earth gravity, we can move the sample cell to 
the bottom end of the solenoid, where the magnetic force is antiparallel to 
the Earth's gravity (see Fig. lb). For even larger effective gravities or studies 
of other fluids, such as water, we may also utilize the new Variable Gravity 



Testbed Facility at the Jet Propulsion Laboratory 15], which is an apparatus 
similar to ours with a higher maximum magnetic field and a room-temperature 
bore. 

Support for this work was provided by the California Institute of Technology's 
Jet Propulsion Laboratory, Low Temperature Science and Quantum Sensor 
Group under contract from the National Aeronautics and Space Administra- 
tion. Support for M.E.T. was also provided by the California Institute of Tech- 
nology Summer Undergraduate Research Fellowships, Tab and Keri Stephens 
Fellowship. 
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Fig. 2. Panel a: a schematic of the experimental setup (not to scale). The sample cell 
containing the liquid oxygen (LOX) consists of a hollow quartz cylinder, stainless 
steel top plate, and a copper bottom plate. A thin foil heater is attached to the 
bottom of the copper plate, outside the sample cell. Resistive temperature sensors 
(red boxes) are located throughout the sample cell with an additional sensor embed- 
ded in the copper plate. Hall effect magnetic field sensors (blue circles) are located 
at either end of the sample cell, along the magnet axis. The cell is located in the 
vacuum space of the cryogenic probe, which is itself immersed in the liquid helium 
used to cool the magnet. Panel b: a block diagram of controls. The pumplines are 
indicated by solid lines. Electronic connections are shown by thick dashed lines. 
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Fig. 3. A typical data run at a saturation pressure of 760 torr and an effective 
gravity of 6.0g. The horizontal axis (common to all three frames) shows elapsed 
time in seconds. The vertical axes show (from top to bottom) applied heater power, 
heater temperature, and system pressure. Some time regions of interest are sepa- 
rated by dashed lines and labeled in the bottom frame: a) warm-up and stabilization, 
b) transient behavior (shaded region), where the onset of boiling causes pressure 
and temperature oscillations until the pressure controller stabilizes, and c) data 
collection. The horizontal dotted line in the center panel indicates the saturation 
temperature. There are no pressure data in region "a" because the valve to the 
sample cell is closed during this phase of the experiment. 
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Reduced superheat t = [T-T b )/T b Reduced superheat t = (T-T b )/T b 

Fig. 4. Heat transport data, extracted from time series data similar to those in Fig. 
3, shown on a linear plot (panel a) and the same data on a log- log plot (panel b). The 
reduced superheat, calculated from the heater temperature T and the saturation 
temperature T&, is displayed on the horizontal scale. The vertical scale shows the 
applied heat flux. The data have been offset vertically in both plots for legibility( 
steps of 1000 in panel a and factors of 10 in panel b). It should be noted that all 
of the fitting curves pass through q = 0, t = if the offset is removed. The symbols 
denote g e g: closed polygons g c g = lg, open polygons g e $ = 6-0g, crosses g e g = lQg. 
The solid lines indicate fits to the data given by Eq. 4 and the parameters in Table 
1. 
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Fig. 5. The sequence of heater powers in these data was chosen to demonstrate the 
memory effect on the heater temperature, as discussed in the text. In panel a, the 
power is alternated between a nominal maximum value of 5 W and a sequence of 
decreasing values. Note that at each return to 5 W, the heater temperature returns 
to the same value of 91.70 ± 0.01 K, approximately shown by the dashed line. In 
panel b, the power sequence is inverted, starting at a nominal minimum power of 
2 W and alternating with an increasing sequence. In contrast to the data in panel 
a, these data do not return to the same heater temperature when the heater power 
returns to 2 W, at the areas indicated by the arrows. 
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